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Abstract: We combine random walks, growth and decay, and convection, in a Monte 
Carlo simulation to model ID interface dynamics with fluctuations. The continuum limit 
corresponds to the deterministic Fisher equation with convection. We find qualitatively 
Q\ • the same type of asymmetry, as well as velocity difference, for interface profiles moving in 
opposite directions. However a transition apparent in the mean-field (continuum) limit is 
not found in the Monte Carlo simulation. 

The Fisher equation [1] combines growth, decay, and diffusion, and was originally 
f^. . used to model population growth subject to limited resources. More recently, a nonlinear 
convective term found in several other pattern-forming systems [2,3] was added to this 
equation [4]. This convection is phenomenologically justified in the case where an external 
l/"") ■ field has a non-zero component parallel to interface motion, and competes with inertial 
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effects. In one dimension, the new equation, which is called FEC (Fisher Equation with 
. Convection), then reads: 
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where all quantities are dimensionless, \i is a positive parameter which cannot be scaled out 



and serves to tune the relative strength of convection, and the density u will be between 



and 1 (thus the convection is said to be towards the right; this is different from advection 
to the right, which is a linear process). The "mean field" equation (1) neglects random 
fluctuations, which we include via a Monte Carlo simulation. Also, the nonlinear nature 
of convection makes its modelling via Monte Carlo non-trivial, an interesting problem in 



itself. Ref. 4 gives more background literature on the Fisher equation and the FEC. 

The Monte Carlo steps are taken as follows. We take Wj to be an integer multiple of 1/n 
(e.g. n=4), so that each site i of our one-dimensional chain carries Uin particles. We take 
into account the diffusion, ^d 2 u/dx 2 , by a random walk. For the growth and saturation, 
u(l — u), we increase Ui with probability 1— uf. each of the Ui n particles produces another 
particle with probability 1 — ui. Thus "saturation" means that it becomes more difficult 
to produce new particles as the ui = 1 value is approached, as when resources are limited. 
A negative probability corresponds to decreasing itj with probability m — 1 to the excess 
over 1 (although we instead set Ui = 1 for speed; this does not change the results). The 
nonlinear convection term, /j,udu/dx, is modelled by moving each particle at site i to the 
right with the density- dependent probability (iUi+i/2, thus introducing a non-linear bias 
in the random walk. Making the probability depend on the right neighbor rather than on 
the site itself eliminates artificial bias due to discretisation. If \x > 1 we rescale the time t 
by a factor \x and divide each of the above three probabilities by \x. Denoting the change in 
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density between two time steps at cell i by Aui, the master equation can then be written 



Am = -(u i+ i - 2m + Ui-i) + Ui(l - Ui) - put — 1 -^~2 — ( 2 ) 

The first term is recognized to be the discretized diffusion, the second one the growth- 
saturation, and the last one is the symmetrized version of the discrete convection (this 
eliminates artificial bias due to discretisation). Note that the effect of convection will be 
lost if n is only 1: it will in effect be simple linear advection. It will also disappear when 
each particle in the system becomes surrounded by empty cells, as is the case when there is 
no growth process (e.g. only diffusion and convection, as in Burgers' Equation [3]). Hence, 
contrary to diffusion, this algorithm will not yield convection if only one particle is in the 
system or the system has only 2 levels (0 and 1, i.e n = 1), hence we choose n > 2. 

One time unit in our one-dimensional simulations (one Monte Carlo step per site) 
corresponds to one convection, one diffusion, and one growth attempt (sequentially), each 
handling all particles in the system at once. The diffusion update was done randomly. The 
convection was tried sequentially right to left, and randomly, without noticeable difference. 
Growth was done sequentially left to right (although, due to the local nature of growth, 
one does not expect this to be important) . We started with a nonzero density in the center 
of the chain and stopped the simulation when the boundaries became occupied. Averages 
were taken over 1000 or 2000 runs. For random numbers we used the built-in 48-bit linear 
congruential function (LCF) on our Hewlett-Packard workstation, and verified consistency 
with the GGL 32-bit LCF, and with a shuffler routine [5] to make the cycle of the random 
number generator effectively infinite. No differences were found. 

Fig.l and Fig. 2 respectively show the density profiles for n = 4 and n = 64, both 
for fj, = 4, at large enough times that the fronts do not change shape. The profile is 
not symmetric, as expected due to the presence of convection, which breaks left-right 
symmetry: the interface thickness on the right is smaller than on the left end of the 
density plateau. However, this asymmetry is less pronounced than in the solution of eq(l) 
in ref. 4. The long-time velocities with which the left and the right interfaces propagate 
are plotted in Fig. 3. This may be compared with the mean-field predictions [6]: \/2 (for 
left front, independent of /i; for right front, for /j, < y/2) and 1/^ + ^/2 (for right front 
when fi > V2). The important qualitative difference is that in the continuum limit the 
two fronts have the same speed \/2 for a certain range of ji < y/2, while in the discrete 
Monte Carlo modelling of convection, the two speeds appear different for all non-zero fi. 
Understanding the reason for this difference will require further investigation. Often mean 
field theories predict transitions that disappear when fluctuations are included (e.g. one 
dimensional Ising model). 

We found that reducing the probabilities for diffusion, growth, saturation and convec- 
tion by the same amount restores the strong asymmetry found in ref.4 between the left and 
the right interface thickness. This is equivalent to decreasing the time step, thus making 
the time evolution of the system less discontinuous. This also allows n to span a greater 
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domain without having to do any time rescalings. The results in Fig. 4 were obtained by 
reducing probabilities by a factor of 1/8 (so \i could go continuously from to 8) and 
looking at \x = 8. The discreteness of the time step then seems crucial in determining 
the amount of asymmetry convection can give to the profile. No qualitative change was 
observed in the /(/-dependence of the asymptotic front velocity, however. Elsewhere [6] we 
have shown that in the corresponding mean field description, the asymptotic velocity is 
constant for \i < 0.5. In our present Monte Carlo simulations, it is not clear whether the 
velocity is really constant or just slowly varying for \i < 0.5. A more thorough investigation 
must be done for small fx. 

To summarize and conclude, fluctuations were introduced in a phenomenological 
model involving diffusion, growth and saturation, and non-linear convection. The con- 
vection seems to have been properly modeled in that basic qualitative features of the 
continuum description are observed both in left-right asymmetry of interface velocities as 
well as interface profile. The fluctuations become less important when a smaller time-scale 
is used (obtained by reducing the probabilities by some factor larger than 1). A more 
thorough evaluation of the effect of discreteness, whether improvements can be made to 
the modelling of convection via Monte Carlo, as well as studies of a Langevin-type equa- 
tion (eq.l with noise) and generalization to higher dimensions, are possible venues worth 
exploring. 
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Fig.l: Density profiles for chains of length 500, with fj, = 4 and n=4 (average over 
1000 realisations). 

Fig.2: As Fig.l, for n=64. 
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Fig. 3: Asymptotic velocities of the left and right interface for n=4 and n=64. Rescaled 
mean field results for eq(l) are shown (right front, solid line; left front, dashed) for quali- 
tative comparison, as they give velocity greater or equal to \f2. 

Fig. 4: As in Fig. 2, but with the probabilities for all four stochastic processes reduced 
by a factor 1/8, and \x = 8. 
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